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ABSTRACT 



The gravitational term for clouds and cores entering in the virial theorem is usu- 
ally assumed to be equal to the gravitational energy, since the contribution to the 
gravitational force from the mass distribution outside the volume of integration is 
assumed to be negligible. Such approximation may not be valid in the presence of 
an important external net potential. In the present work we analyze the effect of an 
external gravitational field on the gravitational budget of a density structure. Our 
cases under analysis are (a) a giant molecular cloud (GMC) with different aspect ra- 
tios embedded within a galactic net potential, including the effects of gravity, shear, 
and inertial forces, and (b) a molecular cloud core embedded within the gravitational 
potential of its parent molecular cloud. 

We find that for roundish GMCs, the tidal tearing due to the shear in the plane of 
the galaxy is compensated by the tidal compression in the z direction. The influence 
of the external effective potential on the total gravitational budget of these clouds is 
relatively small (up to ~ 15 — 25%), although not necessarily negligible. However, for 
more filamentary GMCs, elongated on the plane of the galaxy, the external effective 
potential can be dominant and can even overwhelm self-gravity, regardless of whether 
its main effect on the cloud is to disrupt it or compress it. This may explain the 
presence of some GMCs with few or no signs of massive star formation, such as the 
Taurus or the Maddalena's clouds. 

In the case of dense cores embedded in their parent molecular cloud, we found that 
the gravitational content due to the external field may be more important than the 
gravitational energy of the cores themselves. This effect works in the same direction 
as the gravitational energy, i.e., favoring the collapse of cores. We speculate on the 
implications of these results for star formation models, in particular that apparently 
nearly magnetically critical cores may actually be supercritical due to the effect of the 
external potential. 

Key words: Galaxies: kinematics and dynamics - ISM: general - clouds - kinematics 
and dynamics - Stars: formation 



1 INTRODUCTION 

All known star formation in the Galaxy occurs within dense 
cores in molecular clouds (MCs). The detailed physical na- 
ture of such molecular clouds and their dense cores has been 
a matter of debate over the years. In particular, the su- 
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personic lin ewidths of CO obser ved in molecular clouds for 
first time bv l Wilson et~ai1 (|l970> ) was suggested as in d icativ e 



Goldreich fc Kwanl (|l974h . 
1974) subsequently argued 



of gravit ational contraction by [j 
However, IZuckerman fc Evansl |l 
that, if MCs were collapsing, the star formation efficiency 
should be much larger than observed, and that the molec- 
ular gas in our galaxy should already be exhausted. These 
authors suggested, instead, that the supersonic linewidths 
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were evidence of supersonic turbulence. This idea was widely 
accepted, and turbulence was assumed to be a key ingre- 
dient of molecular cloud sup port against self-gravity. For 
instance, Ide Jong et ail (|l980h calculated hydrostatic mod- 
els of molec ular clo u ds su pported by turbulent (ram) pres- 
sure, while iLarsonl ()l98ll ) found two scaling relations for 
atomic and molecular clouds which were compatible with 
them being gravitationally bound and in approximate virial 
equilibrium. Since then, many observational studies claim 
that molecular clouds and their cores are close to energy 
equipartition (most often referred to as "virial equilibrium" , 
altho ugh this is not necessarily so; see iBallesteros-Paredesl 
2006) between self-gravity, kinetic energy, and, when the 
measurements were avail ab le, the magnetic energ y (e.g., 
iMvers fc Goodmanlll988al lbl; lHeiles fc Troland 12003 ). While 
some authors supported the idea that MCs and t heir cores 
are self-gravitating fe.g. jMcKee and Zweibell992l. although 
with some exceptions, see iBertoldi fc McKed ll 1992)) other s 
argued that MCs may not be self-gravitating ( Blitzl fl994l ) . 
and w ill "need" an ex ternal pressure, in order to be confined 
(e.g., lMalonevlll988h . 

With the development of numerical simulations 
the last decade, it became clear that molecular 



clouds and their cores can exhibit ILarsonl l|l98ll) 



type 



1997; 



relationships (IVzquez-Semadeni et al.l 
iBallesteros-Paredes fc Mac Low! 12002 ) , and nea 
equip artition ( Balles teros-Paredes fc Vazquez-Semademl 
1 1995 ), but they are not in vir ial equili brium 



near energy 



1997 



2004 



( Ballesteros-Parede s fc Vazqu ez-Scma denil 

Shadmehri et al.l |2002| ; iTilley fc Pudrit3 
Ballesteros-Paredes! 120061 ; iDib et all 120071 ). Such clouds 
and their cores can be also more transient than previously 
thought, a fact that is compatible with the lack of post-T 
Tauri s tars (~ 10 Myr old stars) associated to molecular 
clouds l|Bricefio et al. 19971; IBallesteros-Paredes et alj 19991; 



Hartmann et al.l 1200 ll ; IBallesteros-Paredes fc Hartmannl 



20071 ). 



In most virial balance analyses of molecular clouds and 
dense cores, one of the key assumptions is often that the 
gravitational term entering the virial theorem for a cloud or 
core is given only by its gravitational energy, 



£grav = - j p$</\". 



(1) 



where "!> is the gravitational potential produced by the den- 
sity p in the volume V, thus neglecting the gravitational field 
due to the mass distribution outside the volume V. How- 
ever, in principle the mass external to the volume of integra- 
tion can also i nfluence the energy budge t of the cloud/core 
!|Spitzedll97Sl ; lBallesteros-Parededl2006l ). The external dis- 
tribution of mass, together with external forces (such as cen- 
trifugal and Coriolis forces) may either tend to tear apart or 
to compress the cloud/core under analysis. 

The goal of the present work is to numerically eval- 
uate how realistic is the third Viria l Theorem "Myth" dis- 
cussed by IBallesteros-Paredesl (|2006l ) . namely, that the grav- 
itational term can be approximated by the gravitational 



1 By transient we mean clouds that either re-disperse into their 
environment, or collapse to form stars, but that in either case do 
not last much longer than their free-fall times. 



energy. We aim to compare the gravitational energy with 
the effect of external forces for the cases of giant molec- 
ular clouds (GMCs), as well as of the dense cores within 
them. Thus, we evaluate the full external potential term 
in the virial theorem for both (i) a typical GMC within the 
Galaxy, and (ii) a typical dense core within the gravitational 
potential of its parent molecular cloud. 

The organization of the paper is as follows. In fj2]we pro- 
vide a physical interpretation for the external gravitational 
term W ex t, and its contribution to the energy budget of a 
volume of interest. In i|3]we discuss the cases under analysis, 
and in Sjjjwe present our results. Finally, in SjS]we summarize 
our results, and discuss their physical implications, empha- 
sizing their consequences on models of star formation. 



2 TIDAL ENERGY ENTERING THE VT 

The gravitational term entering the VT, given by 
9$ 



IT 



(2) 



is usually taken to be equal to the gravitational energy of 
the cloud, i.e., 

W ~ Sgrav = 1/2 j p $cloud dV, (3) 
J V 

where p is the density, $ is the total gravitational potential, 
and 



'I'cloucl 



-G 



x-x'| 



dV 



(4) 



is the gravitational potential of the cloud, i.e., due only to 
the mass distribution inside the volume V of the cloud. This 
approximation is valid only if the cloud is isolated, i.e., if 
the forces from external agents ar e negligible compared to 
the s elf-gravity of the cloud (e.g., IChandrasekhar fc Fermil 
1953). However, MCs in disk galaxies are confined to 
the midplane [z ~ ±50 pc). Within these disks, most 
of the mass in molecular gas is well organized into spi 



ral arms, bars and/or rings (e.g., I Young fc Scov ille 1991, 
Downes et aT1ll996l: lLoinard et al.l L1996I. Il999l: iDame et all 



200 ll : iHelfer et all 120031 ) . following the gravitational po- 



tential of the old stars, which dominates the dynamics 
ijBinnev fc Tremaind [19871 1 . These patterns suggest, thus, 
that the Galactic gravitational field may be playing a crucial 
role in the global gravitational budget of MCs. An impor- 
tant issue is thus to determine the effect of tidal forces on 
the clouds and their effect on preventing or inducing their 
collapse. 

In order to calculate the contribution of the external 
potential to the energy budget of the cloud, let us separate 
the total gravitational potential into its component due to 
the mass of the cloud alone plus all external agents: 



$ = $cloud + $o 



(5) 



where 3> ex t is the potential produced by all the mass out- 
side the volume of the cloud. Thus, the gravitational term 
entering the VT can be written as: 



W = Bgrav + Wcxt 

where Wext is given by 



(6) 
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Xi p 



dxi 



dV. 



(X) 



This term represents the contribution to a cloud's gravita- 
tional budget from tidal forces. Note that we have defined 
W cx t with a minus sign, such that, if it is negative, it con- 
tributes to the collapse, just as the gravitational energy, but 
if it is positive, it acts against the gravitational energy by 
contributing to the disruption of the cloud. Furthermore, 
its sign is given by the mean concavity of the gravitational 
potential "l? ox t, being positive if the concavity faces down- 
wards, and negative if it faces upwards. To visualize this, 
in Fig. [T] we draw schematically four different situations in 
which a spherical cloud is embedded in an external gravita- 
tional potential <E> ex t. In the first case, (top-left panel) the 
forces (F — — V3>ext) on the right-hand side of the cloud 
are stronger than the forces on the left-hand side. This is 
schematically represented by a darker and larger arrow on 
the right-hand side of the cloud, and with a larger absolute 
value of the slope of the external gravitational potential. 
Since the slope of $ C xt is negative, the forces are positive, 
but stronger at the right- than at the left-side of the cloud. 
The net effect of such an external field is towards disrupting 
the cloud. In fact, for the left-hand side of the cloud, Xi < 0, 
d^cxt/dxi < 0. The product of these quantities is positive, 
and the minus sign of W cx t (eq. [7} gives a negative result for 
Wcxt- However, for the right-hand side of the cloud, Xi > 0, 
and d&ext/dxi < 0. The product is negative, and the mi- 
nus sign in the definition of W eK % gives a positive value for 
Wext- Given the symmetry of the cloud, and since the forces 
are larger on the right-hand side, the total value of H4xt is 
positive, i.e., the external field acts against the gravitational 
energy. In other words, the cloud is being torn apart by the 
external gravitational field. 

For the second case (top-right panel), the dominant con- 
tribution is that from the left-hand side, and the cloud suf- 
fers a net external compression. In this case, a similar anal- 
ysis gives M4xt < 0. Cases 3 and 4 (respectively, bottom left 
and right panels) can be analyzed in a similar way, and Wcxt 
is negative and positive, respectively. In summary, W cyi t con- 
tains the net effect of the tidal forces over the whole volume 
of the cloud, and its sign is defined by the curvature of the 
gravitational potential field. 

Before ending the present section, it is useful to define 



W4xt,z 

We: 



* cxt,y 



x p d&cxt/dx dV, 



y p 9$cxt /dy dV, 



p d<i> cx t/dz dV, 



(8) 



such that Wext = W cxt ,x + W c ^t, v + Wmt,x, in order to un- 
derstand how the three-dimensional potential works in each 
direction. 



3 MODELING THE GRAVITATIONAL FIELD 
AND THE VOLUME OF INTEGRATION 

We evaluate the virial gravitational term in two main cases: 
a giant molecular cloud (GMC) within a disk galaxy, and a 
molecular cloud core at different positions within its parent 
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Figure 1. The tidal energy, Wcxt, is positive if the concavity of 
the external potential <&ext faces downwards, and negative if it 
faces upwards. 



molecular cloud. In every case, following eq. J5), we will 
assume that the total gravitational potential $ is given by 
the addition of two independent potentials: that produced 
by the mass of the object under analysis (GMC or core, 
which will be called $ c ioud), and the external potential (from 
the Galaxy or the GMC), which we will call $ e xt- 



3.1 Giant Molecular clouds within a disk galaxy 

We first consider the effects of a galactic gravitational poten- 
tial field over a giant molecular cloud. We consider a set of 
prolate clouds with constant density n — 50 cm -3 , and ma- 
jor axis of 25 pc. The aspect ratios for each case are: 10:10 
(spherical case), 10:8, 10:6, 10:4, 10:2, and 10:1. We align 
those spheroids in both the x- and {/-directions, in order to 
see how the cloud's energy budget varies as a function of 
the relative orientation of the spheroids and the spiral arms 
within the galaxy. 

The clouds are located within a spiral galaxy with a 
gr avitational potential fie ld given by the model implemented 
by IPichardo et al.l (|2003h . This model uses an axisymmet- 
ric background potential that assembles a bulge, a flat- 
tened disk with a s cale-h eight of 250 pc, as proposed by 
IMivamoto fc Na gai (1975]), and a massive halo extending t o 
a radius of 100 kpc, as proposed bv lAllen fc S antillan ( 199lf ). 
The main adopted parameters are Ro = 8.5 kpc as the Sun's 
galactocentric distance, and Vo(Ro) = 220 km s _1 as the 
circular velocity at the Sun's position. The total mass is 
9x 1O 11 M0, and the local escape velocity is 536 km s _1 . The 
local total mass density is po = 0.15A/q pc -3 . The result- 
ing values for Oort's constants are A — 12.95 km s _1 kpc -1 
and B — —12.93 km s _1 kpc -1 . The full expressions of the 
axisym metric potential can be found in lAllen fc Santillan! 
(| 19911 ). 

In addition to the axisymmetric potential, we have im- 
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plemented a bi-symmetric spiral mod el constructed with 
oblate spheroids as those proposed by ISchmidtl (1 19561) us- 
ing: th e spiral logarithmic locus proposed bv iRoberts et aD 
(1979). The parameters for the spiral arms based on self- 
consis tency studies and a compilation of observational re- 
sults (jPichardo et al.ll200St ) are: the pitch angle is i p = 15.5°; 
the galactocentric radii at which the arms start and end are, 
respectively r s — 3.3 kpc and r e = 12 kpc; the width of 
the spiral arms is 2 kpc. The density fall along the spiral 
arm is exponential with the same scale-length as the disk 
(2.5 kpc). Finally, the ratio of spiral mass to disk mass is 
Ms/Md = 0.0175, which represents a lower limit on the 
range given by self-consistency analysis. 

Every point in the volume of integration will be ana- 
lyzed in a frame of reference located at the center of the 
cloud, which rotates around the Galactic center. This adds 
a centrifugal potential that produces the force 



Or 



Qgr e r 



(9) 



where fio is the angular velocity of the reference frame that 
rotates with the center of the cloud, r is the galactocentric 
distance of the fluid element under consideration, and e r is 
the unit vector in the direction of the galactocentric radius. 
It is also necessary to account for the Coriolis term in our 
modified potential: 

9$cor 



F CO r — 



dr 



= 2Qqv e 



(10) 



where v is the velocity of the point in the volume of integra- 
tion, measured in the frame of reference of the center of the 
cloud, i.e., in the frame that moves with angular frequency 
Oo With this in mind, the effective potential entering in 
eqs. (0 and is given by 



with a — rGMc/2. Following iGieles et al.l (|2006[). we relate 



the mass and size of the GMC by a 
relatiorQ. 



Larson 1 198ll ) type 



M GMC = 540M e 



( ll>c 



^GMC 



(13) 



which, for Mgmc — 10 4 M Q yields rQMC = 4.30 pc. For 
these parameters, the GMCs central density is 4.13 x 
10 3 cm -3 . 

Since GMCs frequently exhibit elongated or filamentary 
structure, we consider two other simple cases of modified 
Plummer spheres (a detailed analysis including a variety 
of more realistic parent cloud structures will be presented 
in a future contribution). In the first place, we studied the 
superposition of two Plummer spheroids: one centered at 
x — and with mass Mgmc, and the second centered at 
x = Stqmc 5 with a mass 4Mgmc, and a radius twice as 
large, according to eq. 1|13[1 . 

Finally, we modify the lPlummerl (|l91ll ) potential to ob- 
tain an elongated structure, in order to mimic a filamentary 
cloud. In this case, the external potential is given by 

GMgmc 



where 

fGMC 



\/f 2 + a 2 ' 



— ) + 1 — 



-, 1/2 



(14) 



(15) 



r x = r G Mc/ 2/3 and r y = r z = r G Mc/~ 1/3 , where / is the 
aspect ratio of the resulting spheroid. For / = 10, r x = 20 pc 
and r y — r z — 2pc. The parameters Mgmc an d »"gmc are 



related also by eq. {13 



V$ cf f = V$gal + V^ccntrif + V$cor 

= V$ ga i - {fl 2 r + 2n v) e r . (11) 

Note that this potential assumes perfect circular orbits, 
while actual orbits are not necessarily perfect circles. How- 
ever, the velocities in both cases are similar within a factor 
not larger than 10%. Since the velocities entering in eq. (|lip 
affect two of the three terms, we do not expect changes larger 
than 10%. 

3.2 Dense cores within a molecular cloud 

In addition to the molecular cloud embedded within a spiral 
galaxy, we analyze the gravitational budget of smaller struc- 
tures: dense cores embedded within its parent GMC. Al- 
though molecular clouds and their cores ex hibit substantial 
subst ructure over a wide range of sizes (e.g., iFalgarone et al.l 
we will consider, for simplicity, that our test core is 
spherical, with constant density (n = 3500 cm" 3 ), and size 
flcore = 0-5 pc, which are typical of dense dark cloud cores 
l|Troland fe Crutcherll200sl ) . We also model the gra vitational 
poten tial field of the parent molecular cloud with a lPlummerl 
j|l91lT) potential, as used bv IGieles et all (I2006D . who ana- 
lyzed the evolution of a stellar cluster interacting with a 
giant molecular cloud. Such potential is written as 

CXt ~ / 2 I 2 ' *• > 

\fr 2 + a z 



4 RESULTS 

4.1 The molecular cloud within a disk galaxy 

Figures [5] and [3] show the ratios W ey ± /F grav through the 
galaxy at z — for our prolate spheroidal GMCs (n = 
50 cm -3 , long axis i2 c loud = 25 pc on the plane of the galaxy, 
and aspect ratios of 10:10, 10:8, 10:6, 10:4, 10:2 and 10:1, 
along the x- and y-directions, respectively). 

The first important result is that the term due to the 
external potential, W cyL t, does not seem to be negligible in 
most of places within the galaxy. This is in clear contrast 
with the typical assumptions in the literature, where the 
gravitational energy -E gra v is thought to be the quantity to 
be compared to other kind of energies in order to explain 
collapse, support, or expansion. Secondly, the maximum and 
minimum possible values of the ratio W e ^t/E SIBX are quite 
similar for both (x- and y- alignment ) cases, although the 
structure of the images differs substantially, specially in the 
cases with large aspect ratios (panels d,e,f in both figures). 

2 Note that this relation has been highly questioned, since 
it seems to be more a consequen ce of o bserva t ional biases , 
rather than a true relation s hip jK egej I l989l: ISc alol Il 990t 
IVzouez-Semadeni et all 1 19971 ; iBallesteros-Paredes fc Mac Low! 
120021 ). However, for the sake of simplicity, the models of clouds 
that we analyzed follow this relationship. 
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Figure 2. Grayscale map of the ratio Wext / -Egrav for GMCs (represented as prolate spheroids) aligned with the x axis. A central circle 
(r ^ 3 kpc) is excluded from the calculation since there arc some missing physical ingredients, such as a galactic bar potential. The 
dashed line represents the W ex t = level. Each case is a prolate spheroid, with largcr-to-smaller ratio of (a) 10:10 , (b) 10:8, (c) 10:6, (d) 
10:4, (e) 10:2 and (f) 10:1, respectively. Note that, for clouds with small axis larger than 10 pc (cases a, b, c, and d), the ratio is positive, 
i.e., W ex t acts in the same sense as gravity. This is due mainly to the strong curvature of the gravitational potential along the galactic 
height z, which dominates in compressing the cloud against any tidal disruption imposed by the shear and/or tidal streams from the 
spiral arms. However, for clouds elongated in the plane of the galaxy by a factor of ~ 10:4 or more (cases e and f), either tidal stretching 
or compression can dominate the energy budget, depending on the position and orientation of the filament in the galaxy. This result 
suggests that the intensity of star formation of a cloud may have a strong dependence on its morphology, its position and its orientation 
within the galaxy. 
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Figure 3. Same as Fig. [2] but for spheroids aligned with the i/-axis. 



We also note that, for the more roundish clouds (cases 
a, b, c, d), the effective potential of the galaxy works in the 
same direction as the self-gravity, giving positive values of 
the term W4xt /-Egrav . This is so because, while the effective 
potential imposes a shear that tries to disrupt the clouds 
in the plane of the galaxy (i.e., W cx t,x + W Q xt, y > 0), the 
gravitational potential of the galaxy along the z axis has a 
confining effect, and it may dominate if the cloud is long 



enough along z compared to its dimensions on the plane of 
the galaxy. The situation is, however, somewhat different 
for more elongated clouds that do not extend very much in 
z (panels d e, f). In these cases, the effective gravitational 
field of the rotating cloud may work either by compress- 
ing (darker zones) or disrupting (light zones) the cloud, de- 
pending on the region where the cloud is located within the 
galaxy. Furthermore, by comparing Figs. [2] and [3] it can 
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be noticed that if the filamentary cloud is perpendicular to 
the radius of the galaxy, the external field works toward 
compressing the cloud (darker zones), while if it is parallel 
to the radius of the galaxy, the external field tries to dis- 
rupt it (light zones), even if it is close to the spiral arms. 
This result is due to the differential rotation, since a cloud 
aligned with the radius of the galaxy will feel the galactic 
differential rotation much more stronger than a cloud that 
is perpendicular to it. 

Finally, it is interesting to note that, for the roundish 
cases, the external contribution to the gravitational budget 
is not larger than ~ 25%. However, for the more elongated 
cases, this contribution can be a large factor, suggesting that 
some filamentary GMCs may actually be gravitationally un- 
bound. This could be the case, for instance, of those few 
GMCs that exhibit no signs of massive star formation, yet 
have masses of the order of 10 4 -10 5 M Q (e.g., Taurus, G 216- 
2.5, Coalsack). 

One may ask how much these results can vary with the 
parameters of the cloud (size and density). Let us assume 
that the cloud is "small enough" , so that the external forces 
over the cloud can be approximated by linear functions of 
the distance to the center of the cloud, then V<I?ext oc -Rcioud- 
In this case, W ox t oc /0-R c i oud , while E gI3lV oc M 2 /R cloud oc 
P 2 -Rcioud- Then, the ratio W ex t/i?grav oc l/p c loud- Thus, for a 
given shape of the cloud, the ratio W EX t / E gra _ v can be scaled 
to different densities as l/p c ioud- In our case, "small enough" 
is smaller than ~ 50 pc. 

In the present work, we have decided to keep the density 
constant and vary the shape of the clouds. However, it is pos- 
sible in principle that the resulting increase of W c ^t/E gla , v 
for spherical clouds for elongated clouds could be due to 
a decrease in the clouds' mass, and therefore their gravi- 
tational energy as the shape becomes more elongated, giv- 
ing the impression of an increase in the relative importance 
of Wext- Therefore, another plausible scenario to explore 
would be to fix the mass, or the column density, while 
varying other parameters of the cloud. However, our re- 
sults can be translated to those cases. For instance, since 
Wext/Sgrav oc 1/pdoud , by keeping constant the mass of the 
clouds, the ratio W ex t/E gra , v grows with size as -Rci oud - In the 
case of constant column density, the Wext/-Egrav oc -R c i ou d. In 
order to show that this is the case, in Fig. [4] displays three 
panels with the ratio W e xt/E eTa , v for spherical clouds as a 
function of x ^ (y — 0, z = 0) for clouds with different 
densities but same size (left panel), different sizes but same 
mass (middle panel), and different sizes, but same column 
density (right panel). It can be seen that W EX t/ E glav oc 1/p 
for all cases (and then oc R 3 and oc R for the constant mass 
and column density respectively). Thus, for a given cloud, 
one can calculate its mass, or its column density, and then, 
calculate how much the ratio Wext / -E grav will vary if the size 
varies too. 

We finally note that our results do not change strongly 
by placing the clouds above or below the plane of the galaxy. 
At first glance, one might think that a cloud located 50 or 
100 pc away from the midplane will experience a tidal stress 
along z, rather than a tidal compression. However, we note 
that the gravitational potential along z has always positive 
concavity, and then the situation is just like the one shown 
in Fig[T] cases 2 or 3. It is important to note, however, that 
the radius of curvature of the gravitational potential grows 



with z, and thus the vertical compression decreases as z 
increases. 

4.2 The core within a GMC 

Figure [5] shows both the potential (upper panel) and the 
ratios W c *t/E grav , W BXt , x /E gla , v , and W c *t, y /E grav (lower 
panel) as a function of distance to the center of the cloud, 
for our dense core (n — 3500 cm -3 , i? C orc = 0.5 pc) em- 
bedded in a single Plummer potential. As discussed in 321 
the curvature of the potential defines whether W is posi- 
tive or negative. For instance, W cx t,x is negative whenever 
$ C xt has upwards concavity, and positive wherever <E> e xt has 
downward concavity. On the other hand, W ex t, y is always 
negative at every position along the x axis because, along 
the radius, $ ex t has upwards concavity in the azimuthal di- 
rection. 

More important, however, is to note in this figure that 
W ox t grows substantially when the core is located in the 
inner parts of the molecular cloud. In fact, it is larger (in 
absoluve value) than -E gra v at the central parts of the poten- 
tial well, enhancing the importance of accounting for Wext 
when analyzing the energy budget of dense cores. 

As a step towards taking a more realistic potential for 
a molecular cloud, we show in Fig. the potential (up- 
per panel) and the ratios Wcxt/Egrav, W c ^t,x/E gI!3lV , and 
Wext.^/Sgrav (lower panel) for the same dense core embed- 
ded in the potential produced by two Plummer spheroids, 
centered at different locations (x = y — zO pc), and 
x — 12.5, y — z — pc), as a function of x. In Fig. [7] 
we show a grayscale x — y map of the ratio Wext/-E grav on 
the (x, y) plane. As in the previous case, the external field 
again contributes substantially to the energy balance of the 
core, especially near the minimum of the potential well of 
the cloud. It should be noticed also that what defines the 
relative importance of W cx t is the curvature of the potential 
field, i.e., the width of the potential, and not only its depth 
(or height). For instance, the Plummer sphere centered at 
x = 12.5, even though it is more massive, is also more ex- 
tended, producing a potential well with a larger radius of 
curvature. This translates into a lower value of W ex t at the 
position of the central region of that sphere, compared to 
the contribution to W ex t of the other Plummer sphere. 

One may ask, in general terms, if this potential has 
places where it can produce disruption or not. To answer 
this question, in Fig. [5] we show the potential field of the last 
example (two Plummer spheroids) , both in grayscale and in 
a surface representation, as a function of x and y. First of 
all, note that the x axis spans up to 20 units, while the y 
axis spans only 8. Second, note that the plot is symmetric 
with respect to the axes x = and y — 0. Third, note also 
that the dependence on z (not shown here) is the same as 
the dependence on y. 

According to our discussion above, disruption occurs in 
regions with predominant downwards concavity. In our ex- 
ample (see Fig.[T]), the regions that are potentially disruptive 
are (a) close to x = 4, y — 0, z — pc, and (b)regions at large 
x, y and z. However, as can be seen in Fig. [71 W cx t is always 
compressive (Wext / '-Egrav > everywhere). The reasons for 
this are that (a) At large y (or z), the curvature along y is 
positive (upward concavity). However, the curvature along 
x is negative, and with a smaller radius of curvature (see 
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Figure 4. Variation of the ratio W ex t/ E grBjV for spherical clouds of the same size but different densities (left panel), same mass but 
different sizes (middle panel), and same column density but different sizes (right panel). Note that W e nt/E gT!lv °c 1/p, for all cases, which 
implies that Wext / -E gra v oc R^ loud and Woxt / E grav oc R c i ou d f° r the cases of constant mass and column density, respectively. 



F 




-20 t 
1 .5 F 




-0.5 r , , , . . . . . . . . . . . J 
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x [pc] 

Figure 5. Upper panel: Plummer's potential as a func- 
tion of the radius, which we take along the x axis. Lower 
panel: Wext/Bgrav (solid line), W ex t,x/E gT!lv (dashed line) and 
Wext,t//£-grav (dotted-dashed line) ratios for a R CO re = 0.5 pc, 
n = 3500 cm" 3 core embed ded in a molecular cloud of f 4 Mq 
following a lPlummerl &1 91 ll) potential. Note that the tidal contri- 
bution to the energetic content of the core (Wcxt) are as important 
as the self-gravitational energy (E gIav ) of the core. 

Fig. [8] keeping in mind the different ranges spanned by the 
axes). Thus, the net effect of V^cxt over the cloud is com- 
pressive, (b) The local maximum at x ~ 4, y = 0, z = is 
actually a saddle point, and the concavity along x, although 
negative, has a larger radius of curvature than the concavity 
along y, which is positive (again, note that the y axis is more 
compressed than the x axis). Thus, in both (a) and (b) cases, 
the positive curvature dominates, producing a positive value 
for Woxt/Ugrav Then, this potential produces compressions 
everywhere, at least for spherical cores. 



I 
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Figure 6. Same as fig. [5] but for a dense core embedded in a 
gravitational potential given by the superposition of two Plummer 
spheroids. Dashed lines in the upper panel show the potential of 
the individual spheroids. 



Let us now consider an elongated molecular cloud, for 
which the potential is given by a highly elongated Plummer 
spheroid that has an aspect ratio of 10 f ^3.2[) . Fig. [9] shows 
the ratio W cx t / E SI!iv in this case. The dotted line marks the 
place where the ratio is equal to zero. We note that the 
maximum contribution from Woxt to the budget is now sub- 
stantially larger (~ 3.5 times the value of the gravitational 
energy -E gr av). This is due to the fact that the density of 
a filament changes on short scales along the direction per- 
pendicular to its long axis, producing a potential well with 
small radius of curvature in this direction. Note also that the 
upper-left region, (above the dotted VK ex t/i?grav = line), 
has negative values of the ratio Wext /-Egrav , i.e., positive val- 
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Figure 7. Grayscale two-dimensional version of the lower panel 
of Fig. \E\ Note that the scales in y are expanded by a factor of 
2.5. The ratio Wext/Eg has cylindrical symmetry around the x 
axis. 




Figure 8. Gravitational potential produced by the superposi- 
tion of the two Plummer spheroids in the x — y plane, both as 
a grayscale and surface representation. Both minima occur at 
y = 0,z = 0, while at x ~ 4 pc occurs a saddle point. A core 
located in this point will have strong compression if it is long 
enough along y, compared to its x-size, but strong disruption if 
it is long enough in x compared to its y-size (note that the x axis 
spans up to 20 units, while the y axis spans only 8.) 



10 20 30 

x [pc] 

Figure 9. Grayscale of the ratio Wcxt/^grav for a Plummer elon- 
gated spheroid, with aspect ratio 10. The dotted line marks the 
Wcxt = level (notice that the y axis is expanded by a factor 
of 10). Elongated structures seem to produce larger values of the 
ratio W e xt/E 

grav • 

ues of Woxt. Although those values are small, the existence 
of this region shows that there could be places where an iso- 
lated core, in principle, might be disrupted by tidal forces 
from a larger cloud. However, their small absolute values 
imply that, if such a core is in energy equipartition between 
self-gravity and other forms of energy (e.g., magnetic, turbu- 
lent, thermal), the tidal disruption will play a negligible role, 
unless the core is substantially larger than the one modeled 
here. 

In all the cases shown in the present section, the vari- 
ation of the gravitational potential is soft enough that the 
external force V3> C xt can be considered linear through the 
size of the core. Thus, as in the case of the GMC embed- 
ded in a galactic gravitational field, W ox t/E gla , v oc l/p, and 
our results can be easily translated to cores with different 
densities. The case of a core within a more realistic parent 
cloud (e.g., a filamentary GMC from numerical simulations) 
is deferred to a further contribution. 

4.3 General remarks 

In summary, for GMCs extended large enough in z within 
a galactic potential, the external term Woxt is negative al- 
most everywhere because the compression along z domi- 
nates. However, for clouds confined to the plane of the galaxy 
(for sizes smaller than 10 pc in z), the galaxy may either 
compress or disrupt them, depending on the relative orien- 
tation between the cloud and the direction of rotation of 
the galaxy. The spiral arms tend to produce compression, 
although this also depends on the orientation and the size 
of the elongated cloud. 

In the case of the cores, where we have only consid- 
ered the effects of the gravitational potential of their parent 



10 Ballesteros-Paredes, et al. 



cloud, the tidal energy W G xt of the structure embedded in a 
bigger potential can be important close to the local minima, 
and acts in the same sense than the gravitational energy: 
by compressing the cloud. We have not found any case in 
which Woxt is positive and comparable to the gravitational 
energy, suggesting that the tidal energy of a GMC over its 
cores helps in their confinement and collapse, but not in their 
disruption. 



5 DISCUSSION AND CONCLUSIONS 

Ranging from GMC sizes (from a few to several tens of 
parsecs), down to prestellar core sizes 0.1 pc), molec- 
ular clouds and their substr ucture are often thought to be 
in virial equilibrium (e.g., iMcKee fc Ostrikerl 12007. and 
references therein). This state is inferred from a compari- 
son between their internal (kinetic, thermal, magnetic) and 
gravitational energies. In previous work, however, we have 
emphasized the importance of accounting for the role of 
some additional energy terms that should be considered, 
but that a re usually neglected when the virial theorem is 
used (see Ballesteros-Paredes & Vazqucz-Scmadcni 1997; 
Ballesteros-Paredes etal ]|l999a: Ballestero s-Paredes! |200~ 



Dib et al.ll2007l ). In particular. iBallesteros-ParedesI ( 20Qi 
argued that there are at least six assumptions usually made 
in virial theorem analyses that are frequently not fulfilled 
in the interstellar medium. In the present work we have fo- 
cused on the role of the net external potential has on a 
density structure (a cloud or a core). As a first attempt to 
evaluate how realistic this assumption is, we have analyzed 
the contribution of the mass external to the volume of in- 
tegration to the gravitational term W4xt, and compared it 
to the gravitational energy, in two main cases: a GMC em- 
bedded within a spiral galaxy, and a molecular cloud core 
within its parent molecular cloud. 

Although our analysis does not exhaust the possibili- 
ties of cloud and core shapes and configurations, the results 
presented here suggest that by neglecting the influence of 
the net potential external to a given core, observational es- 
timates can lead to wrong estimations of its actual gravita- 
tional content. 

We have found that the gravitational field of a Galaxy 
may have an important influence on the energy budget of 
GMCs. While the galactic rotation tries to tear apart clouds 
that are elongated in the radial direction, the field along z 
tend to compress the clouds. The detailed gravitational con- 
tent of GMCs depends, thus, on their size, shape, position 
and orientation within the galaxy. This could be the case, 
for instance, of Maddalena's cloud, G216-2.5, which exhibits 
no si gns of massive star formation(see, e.g., Mcgcath et al. 
2009, for the current status of star formation going on in 
G216-2.5), in spite of its having a mass ~ 6.6 x 10 s M e . 

The situation for dense cores within simplified poten- 
tials for molecular clouds is similar. The gravitational field 
may tend to disrupt the core in one direction, while com- 
pressing it in another. The contribution of Wext depends 
strongly on the position of the core inside the cloud, and on 
the details of the gravitational potential of the later. In all 
the cases shown here, however, this contribution works in 
the same direction than the self-gravity, i.e., towards com- 
pressing the cloud. Moreover, it can be of the order of mag- 



nitude or even larger than the gravitational energy, specially 
in elongated, filamentary clouds. This is due to the fact that 
filaments can have strong upwards concavity in the direction 
perpendicular to their main axis. 

An interesting result is that tidal disruption tends to 
be more important than tidal compression for larger objects, 
while smaller objects tend to suffer larger compressions. The 
reason for this is that there are no sources of gravitational 
repulsion, only gravitational attraction. Then, the radii of 
curvature when the concavity is positive are smaller than 
when the concavity is negative. This can be pictured by 
imagining the gravitational field as a tense sheet supporting 
weights placed on it. On such arrangement, positive (up- 
wards) concavities at the position of the weights have small 
radii of curvature, while negative (downwards) concavities 
always extend over larger distances. Thus, for a cloud of 
moderate size located where the external gravitational field 
exhibits negative concavity, W ex _t will be positive, but rel- 
atively small with respect to other forms of energy (self- 
gravity, magnetic, turbulent, etc.). However, if such a cloud 
is located in a place where the curvature of the external field 
is positive, its W ex t will be important, since the field varies 
strongly over smaller distances. 

As a conclusion, the results shown here suggest that, 
on large scales, i.e., GMCs within galaxies, the former can 
be farther from real equipartition, either disrupted or com- 
pressed, depending on their shape and orientation within 
the galaxy. On smaller scales, it seems that it is important 
to determine observationally the role of the net external po- 
tential on local dense cores. 
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